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We analyze the chargino contributions to the K — K mixing and e' in the mass 
, insertion approximation and derive the corresponding bounds on the mass 

insertion parameters. We find that the chargino contributions can significantly 
enlarge the regions of the parameter space where CP violation can be fully 
super symmetric. In principle, the observed values of e and e' may be entirely 



Abstract 



due to the chargino - up-squark loops. 



A convenient way to parameterize SUSY contributions to the flavor changing processes 
is to employ the so called mass insertion approximation [Ij]. The advantage of this ap- 
proach is that it allows to treat such contributions in a model independent way without 
resorting to specific assumptions about the SUSY flavor structures (a technical definition 
^ ! of this approximation will be given below). 

The gluino contributions to the kaon observables in the mass insertion approximation 
have been studied in detail @-[|J, but the chargino contributions have not received sim- 
ilar attention. The latter have been considered either in the context of minimal flavor 
violation, that is in SUSY models with the flavor mixing given by the CKM matrix [Q, 
H or as contributing to the K — K mixing only j7|. In general, the flavor structure in 
the squark sector may be very complicated. In particular, flavor patterns in the up and 
down sectors can be entirely different, which may result in the dominance of the chargino 
contributions to the K and B observables. On the other hand, the neutralino contribu- 
tions involve the same mass insertions as the gluino ones (i.e. down type mass insertions) 
and thus cannot qualitatively change the picture. In this Letter, we study the chargino 
contributions to the K — K mixing and e' using the mass insertion approximation and 
derive the corresponding bounds on the mass insertion parameters. 

Let us first consider the K — K mixing. The two observables of primary interest are 
the Kl-Ks mass difference and indirect CP violation in K — > 7r7r decays: 

AM K = M Kl - M Ks , 
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Figure 1: Leading chargino - up-squark contribution to K — K mixing. 
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The experimental values for these parameters are AMjf ~ 3.489 x 1CT 15 GeV and e ~ 
2.28 x 10~ 3 . The Standard Model predictions for them lie in the ballpark of the measured 
values, however a precise prediction cannot be made due to the hadronic and CKM 
uncertainties. 

Generally, AMk and e can be calculated via 

AM K = 2Re{K°\H* s=2 \K } , 

1 -lm(K°\H^ s=2 \K ) . (2) 
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Here H^f =2 is the effective Hamiltonian for the AS* = 2 transition. It can be expressed 
via the Operator Product Expansion (OPE) as 

#e A / =2 = £a(//)a, (3) 

i 

where Cj(/i) are the Wilson coefficients and Qi are the relevant local operators. The main 
uncertainty in this calculation arises from the matrix elements of Qi, whereas the Wilson 
coefficients can be reliably calculated at high energies and evolved down to low energies 
via the Renormalization Group (RG) running. 

In supersymmetric extensions of the SM, the dominant chargino contribution to H^ s=2 
comes from the "super-box" diagram in Fig.l. We perform our calculations in the super 
CKM basis, i.e. the basis in which the gluino-quark-squark vertices are flavor-diagonal. In 
this basis, the chargino - left quark - left squark vertices involve the usual CKM matrix: 

A£ = -9EE V kl K* ba dfio 2 {xt L yu b L , (4) 

k a,b 
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where K is the CKM matrix, a, b are the flavor indices, k — 1, 2 labels the chargino mass 
eigenstates, and V, U are the chargino mixing matrices defined by 



M x 



M-> \f2M w sin (3 ' 

V2M W cos fi 
U* M x+ V- 1 = diag(m x +,m x +) . (5) 



Only the gaugino components of the charginos lead to significant contributions to the 
K — K mixing since the higgsino couplings are suppressed by the quark masses (except for 
the stop coupling) and are not important even at large tan (3. The stop loop contribution 
is suppressed by the CKM mixing at the vertices: each vertex involving the stop is 
suppressed by A 2 or A 3 with A being the Cabibbo mixing, whereas we will be working in 
0(\) order. The super-box involving higgsino interactions with the stops depends on the 
left-right mass insertions and, as will be clear later, does not lead to useful constraints on 
the SUSY flavor structures. 

Due to the gaugino dominance, chargino-squark loops will generate a significant con- 
tribution to only one operator 

q 1 = gfrfdz 4iX , (6) 

similarly to the Standard Model (a, (3 are the color indices). The corresponding Wilson 
coefficient is given by the sum of the Standard Model and the chargino contributions: 

C 1 ( /U )=C 1 (/i) 5M + C 1 (/i)* + . (7) 

Generally, there are additional contributions from gluinos and the Higgs sector, but they 
are not correlated with the chargino contributions and are unimportant for the present 
study. We calculate Ci(/i)* + using the mass insertion approximation. That is, we express 
the left-left squark propagator as 



i 5 ab i (8m 2 ) 



(Hint) = i - m 2 l - SmTJ - + ( ;ri'^2 > (8) 

where 1 is the unit matrix and m is the average up-squark mass. The SUSY contributions 
are parameterized in terms of the dimensionless parameters (5'l L ) ab = (5m 2 ) ab /m 2 . The 
corresponding Wilson coefficient is calculated to be 



where X; L = m 2 + / m 2 and 



2 + 5x — x 2 6x In x 
[l-xf + (l-x) 4 



It is interesting to note that "flavor- conserving" mass insertions (S^ L )aa contribute to 
Ci(Mw), unlike for the gluino case. Such mass insertions arise from non- degeneracy 
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M 2 \ m 


300 


500 


700 


900 


150 


0.04 


0.06 


0.08 


0.09 


250 


0.07 


0.08 


0.09 


0.11 


350 


0.09 


0.10 


0.11 


0.12 


450 


0.12 


0.12 


0.13 


0.14 



Table 1: Bounds on ^Re [(^ L ) 2 i] 2 | from AM K (assuming a zero CKM phase). To obtain 
the corresponding bounds on 8, these entries are to be multiplied by 4.6. These bounds 
are largely insensitive to tan (3 in the range 3-40 and to \x in the range 200 — 500 GeV. 

of the squark masses and are proportional to the difference of the average squark mass 
squared and the diagonal matrix elements of the squark mass matrix. If the diagonal 
elements are equal, the "flavor-conserving" mass insertions drop out of the sum due to 
the GIM cancellations. 

The flavor structure appearing in Eq|| can be expanded in powers of A: 

52K 2 (S2 L ) ab K bl = {5 u LL )n + A [(^x)n - (S U LL )22] + 0(X 2 ) . (11) 

a,b 

Assuming the presence of one type of the mass insertions at a time in Eq.|Tl| at each order 
in A, one can derive constraints on (5^)21 and 8 = (<^ L )n — (8^ L )22 imposed by AMk 
and e. A much weaker constraint on (5^l)3i can a l so t> e obtained if we are to keep 0(X 2 ) 
terms in Eq.|Tl]. 

To derive constraints on the mass insertions, one has to take into account the RG 
evolution of the Wilson coefficients. In our numerical numerical analysis, we use the NLO 
QCD result Ci(//)* + ~ 0.8 C^Mw)^ with fi = 2 GeV §. The matrix element of Q 1 
is computed via (K \Qx\K°) = ^M k Jk b M with the lattice value B^/i) = 0.61 ||. In 
addition, the SM contribution should be taken into account. Its detailed discussion can 
be found in Ref.||. In our analysis, we assume a zero CKM phase which corresponds 
to a conservative bound on the mass insertion. The Wolfenstein parameters are set to 
A = 0.847 and p = 0.4. The other relevant constants are M K = 0.498 GeV and fx = 0.16 
GeV. 

The resulting bounds on (8^)21 an d 8 as functions of M 2 and the average squark mass 
m are presented in Tables 1 and 2. We find that these bounds are largely insensitive to 
tan p in the range 3-40 and to \x in the range 200 — 500 GeV. This can be understood since 
these parameters do not significantly affect the gaugino components of the charginos and 
their couplings. Note that 8 is real due to the hermiticity of the squark mass matrix and 
therefore does not contribute to e. The presented bounds on the real part of {8^)21 are 
a bit stronger than those derived from the gluino contribution to the D — D mixing [[3], 
whereas the imaginary part of (8f jL )2i is not constrained by any other FCNC processes. 

In principle, it is possible to constrain the (5^)31 mass insertion as well. At order 
A 2 , there are two contributions in Eq.|TT|: from (5xz)3i an d (^ll)i2 + (^ll)21- Assuming 
no cancellations between these two terms, the constraints on (8^1)31 are obtained by 
multiplying the bounds in Tables 1 and 2 by (AX 2 )^ 1 ~ 24. Clearly, this leaves the real 
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7.2 x 10~ 3 


9.1 x 10" 


-3 


1.1 x 10" 


-2 


250 


7.8 x 10" 


-3 


9.2 x 10~ 3 


1.1 x 10- 


-2 


1.3 x 10- 


-2 


350 


1.1 x 10- 


-2 


1.2 x 10~ 2 


1.3 x 10- 


-2 


1.5 x 10- 


-2 


450 


1.5 x 10- 


-2 


1.5 x 10~ 2 


1.6 x 10- 


-2 


1.7 x 10- 


-2 



Table 2: Bounds on y |Im [(5^)21] I from e. These bounds are largely insensitive to tan/5 
in the range 3-40 and to fi in the range 200 — 500 GeV. 



part of [(<5i L )3i] 2 essentially unconstrained, while the bound on y Im [(5^)3i] 2 is of order 
10 _1 . We note that a similar constraint on (<^r)i3 can be derived from the higgsino-stop 
contribution, however such a constraint is typically satisfied automatically (especially if 
the squarks are heavy) since (Sl R )i3 ~ em t /m with e < 1 being the 1-3 mixing the 
left-right sector. 

Next let us consider the chargino contribution to e' using the same approximations. 
The e' parameter is a measure of direct CP violation in K — > nn decays given by 



e' 
s 



CO 



V2\e\RjeA 



- (imAo - - lmA 2 ) , 
o V u J 



(12) 



where A ^ are the amplitudes for the AI = 1/2,3/2 transitions and uo = R,eA 2 /ReA ~ 
1/22. Experimentally it has been found to be Re(e'/e) ~ 1.9 x 10~ 3 which provides firm 
evidence for the existence of direct CP violation. This value can be accommodated in the 
Standard Model although the theoretical prediction involves large uncertainties. 
The effective Hamiltonian for the AS = 1 transition is given by 



ttA5=1 



£c^)4 



(13) 



Just as in the Standard Model, two operators, Qq and Q$, play the dominant role. They 
originate from the gluon and electroweak pengiun diagrams (Fig. 2) and are defined by 



Qs = (s a d%„ A ]T (<?q a ) v+ A, 

q=u,d,s 

Qs = l(s a d%„ A £ e q (qPq a ) v+A , 

q=u,d,s 



(14) 



with (ff)v-A = /7/*(l — 75)/- Their matrix elements are enhanced by (m^/m s ) 2 com- 
pared to those of the other operators: 



((7T7T) I=0 \Q 6 \K ) = 

((nn) I=2 \Q 8 \K°) = V3 




m K 



1 2 



m s (fjt) +m d (fj) 
2 



m 2 K (f K ~ U) B 6 , 



m K 



m s (n) +m d (jj) 



m 2 K U B 8 , 



(15) 
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where Bq$ are the bag parameters. In addition, the contributions of these operators are 
enhanced by the QCD corrections. Although the Wilson coefficient of Q 8 is suppressed by 
a/a s compared to that of Qq, its contribution to e' is enhanced by 1/u and is significant. 
In fact, it provides the dominant contribution in our analysis. 




q q q q 



Figure 2: Leading chargino - up-squark contributions to e' (a "mirror" diagram is not 
shown) . 




The relevant QCD corrections in the context of the MSSM with minimal flavor vio- 
lation have been studied in Ref.0 and later, in more detail, in Ref.||. To account for 
a general flavor structure in the mass insertion approximation, only the loop functions 
of Ref. are to be modified. In our numerical analysis we use the parameterization of 
Ref. U and express the chargino contribution to e' / e as 

= Im(^£ K* a2 (6 u LL ) ab K b )j .F £f , (16) 

where 

F £ , = (P x + P Y + P z ) F z + -P z F 7 + P E F g . (17) 

Here we have omitted the box diagram contributions which are negligible ||. The param- 
eters Pi include the relevant matrix elements and NLO QCD corrections, and are given 
by Px = 0.58, Py = 0.48, Pz = —7.67, and Pe = —0.82. The quantities Fi are functions 
of supersymmetric parameters resulting from the gluon, photon, and Z penguin diagrams 
(Fig. 2) and are calculated in the mass insertion approximation. Explicitly, 
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= 2^fEl^l 2 A(x i ) 
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M 2 \ m 
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500 
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900 


150 


0.11 


0.11 


0.13 


0.16 


250 


0.17 




0.87 


0.64 


350 


0.12 


0.29 


0.74 




450 


0.12 


0.23 


0.42 


0.79 



Table 3: Bounds on |Im(52 L ) 2 i| from e'. For some parameter values the mass insertions 
are unconstrained due to the cancellations of different contributions to e'. These bounds 
are largely insensitive to tan (3 in the range 3-40; [i is set to 200 GeV. 
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F Z = --2j>a| 2 ^ 

8 i 



(18) 



where Xi = m 2 + /m 2 and the loop functions are given by 

Ate) 
f { z'(x,y) 



1 - 6x + 18s 2 - 10a; 3 - 3s 4 + 12a; 3 In a; 

18(x- l) 5 ' 
22 - 60a; + 45a; 2 - 4a; 3 - 3a; 4 + 3(3 - 9a; 2 + 4a; 3 ) In x 
27(a; - l) 5 ' 
(J) N (y — 1) [(x — l)(x 2 — x 2 y + xy 2 — y 2 ) + x 2 (y — 1) In a;] — (x — l) 2 y 2 lny 



xy 



16(x-l) 2 (y-l) 2 (y-x) 
(y — 1) [(x — 1) (x — y) + x(y — 1) In x] — (x — l) 2 y\ny 



l(x-l) 2 (y-l) 2 (y-x) 



(19) 



As noted in ref. 0, the dominant contribution typically comes from the Z-penguin dia- 
gram, especially if the SUSY particles are heavy. This can be seen as follows. Due to 
the gauge invariance, the gsidi and jSLdi vertices are proportional to the second power 
of the momentum transfer, i.e. (q^q u — g '^q 2 ) / 'rn 2 '. This momentum dependence is can- 
celled by the gluon (photon) propagator which leads to the suppression factor 1/m 2 in 
the final result. On the other hand, the Zsidi vertex exists at q 2 = due to the weak 
current non- conservation and is momentum-independent to leading order. It is given by 
a dimensionless function of the ratios of the SUSY particles' masses. The Z propagator 
then leads to the suppression factor 1/Mf which is much milder than 1/m 2 appearing in 
the gluon and photon contributions. 

The resulting bounds on Im(5£ L )2i are presented in Table 3. Note that there is no 
SM contribution to e' since we assume a vanishing CKM phase. In the limit of heavy 
superpartners, these bounds become insensitive to the SUSY mass scale. This occurs due 
to the dominance of the Z penguin contribution. Indeed, the contributions of the photon 
and gluon penguins fall off as 1/(SUSY scale) 2 as can be seen from Eq.[L8]. On the other 
hand, the Z-penguin contribution stays constant in the "decoupling" limit. This may seem 
to conflict with the intuitive expectation of the decoupling of heavy particles. However, 
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the proper decoupling behaviour is obtained when the flavor violating parameters 5m 2 ab 
are kept constant (or if they grow slower than the masses of the superpartners) . To put 
it in a slightly different way, the decoupling should be expected when the mass splittings 
among the squarks grow slower than the masses themselves. 

It is noteworthy that (for universal GUT scale gaugino masses of around 200 GeV) 
the bounds on Im(5^ L )2i are slightly stronger than those on Im(5^ L )2i derived from the 
gluino contribution to e' ||. The suppression due to the weaker coupling is compensated 
by a larger loop function mainly due to the presence of the diagram on the right in Fig. 2. 

These results show that to have a chargino-induced e' would require a relatively large 
LL mass insertion (O(10 -1 )) which typically violates the constraints from AMk and e. 
Yet, it is possible to saturate e and e' with the chargino contributions in corners of the 
parameter space. For instance, taking M 2 = 450 GeV and m = 300 GeV, e requires 

2 lm(5 u LL ) 21 Re(5 u LL ) 21 ~ 2.3 x 10~ 4 . (20) 

Then, assuming Im(<^ L ) 21 ~ 0.12 to produce e', Re(5£ L ) 2 i has to be 9 x 10~ 4 . These 
values are in marginal agreement with the AM^ bound: 



[Re(^L)2if-[Im(^ L ) 21 ] 2 | <0.12. (21) 

The main lesson, however, is that combining the chargino and the gluino contributions 
can provide fully supersymmetric e and e' in considerable regions of the parameter spaceQ. 
For example, only a small (Im(5^ i? ) 2 i ~ 10~ 5 ) mass insertion in the down-sector is required 
to generate the observed value of e' ||. Then e can be entirely due to the mass insertions 
in the up-sector: Im(<5^ L ) 2 i ~ 10 -3 and Re(5^ L ) 2 i ~ 10~ 2 . Generally, this does not require 
large SUSY CP-phases and may be accommodated in the framework of approximate CP 
symmetry |]TTJ, which is motivated by the strong EDM bounds (for a review see JT4]]). 



Alternatively, the CP-phases can be order one but enter only flavor-off-diagonal quantities 
which occurs in models with hermitian flavor structures ]T5| . Clearly, the regions of the 
parameter space where CP violation can be fully supersymmetric are significantly enlarged 
if both the gluino and the chargino contributions are included. Of course, it remains a 
challenge to build a realistic well-motivated model with all of the required features. 
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